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Current, power, and attenuation mea- 
surements with rf SQUIDs are based on 
the fact that the voltage from the mi- 
crowave readout circuit can be made a 
nearly sinusoidal function of the mag- 
netic flux threading the SQUID. We 
point out here that an asymmetric dc 
SQUID with sufficiently low inductance 
can achieve a very sinusoidal output 
voltage with good modulation depth. 
The spectral purity of the sinusoid can 
be substantially better than that obtained 
with rf SQUID systems. The purity im- 
proves with increasing asymmetry of 
the junction critical currents, and 



decreasing values of the LI C product, 
where L is the SQUID inductance and 
I c is the smaller of the critical currents. 
Results of several calculations are pre- 
sented. Substantial improvement in 
SQUID methods of rf current, power, 
or attenuation measurement may thus be 
possible with use of such dc SQUIDs. 
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Introduction 

In the early 1970s, Kamper et al. [1-5] 1 recog- 
nized that superconducting quantum interference 
devices (SQUIDs), with their periodic response to 
magnetic flux <J>, had the potential for being used in 
a completely new kind of current, power, or atten- 
uation measurement. In the systems developed, the 
rf impedance of a single-junction SQUID (rf 
SQUID) was measured. The impedance variations 
were reflected into a resonant circuit and measured 
as variations in phase or amplitude of the output rf 
voltage. If the voltage response is purely sinusoidal 
in flux with period equal to the flux quantum 
<& =h/2e, 
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F=F cos(277-<J>/<I>o), 



(1) 



the time-averaged voltage response to a sinu- 
soidally varying flux of amplitude 4> m takes the 
form 



V= V Jo (27T<i> m /<I>o). 



(2) 



Numbers in brackets indicate literature references. 



Equation (1) is achieved approximately with use of 
an overdamped SQUID and a broadband mi- 
crowave readout circuit [4]. By determining the ze- 
roes of the response experimentally as the flux 
amplitude is varied, and equating them to the ze- 
roes of the Bessel function J , one can determine 
<J> m or its equivalent in terms of current. 

The periodic response to an impressed flux is, 
however, not ideally sinusoidal, and this necessi- 
tates an elaborate series of corrective procedures 
[2-7]. The purpose of this note is to point out that a 
dc SQUID (a superconducting loop containing 
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two Josephson junctions — see fig. 1) with asym- 
metric critical currents and sufficiently small in- 
ductance and capacitance can have a considerably 
purer sinusoidal voltage output with adequate 
modulation. Thus one of the major sources of diffi- 
culty in realizing a Josephson junction device for rf 
measurements can potentially be removed. Al- 
though the SQUID inductance is small, the current 
in the input coil necessary to couple flux into the 
SQUID can be kept at a reasonable level. The 
noise in such SQUIDs is also expected to be small. 
The threshold characteristic (i.e., maximum zero- 
voltage current vs. flux) of a symmetric dc SQUID 
is not at all sinusoidal. Figure 2 shows an example 
for a 4-pH SQUID with equal critical currents of 
55 jxA. It has long been known, however, that the 



voltage developed across a symmetric SQUID has a 
sinusoidal appearance when the bias current is suf- 
ficiently greater than the maximum critical current 
of the device. This is observed experimentally and 
from calculations (see, e.g., [8] and references 
therein). 

That an asymmetric dc SQUID can exhibit a si- 
nusoidal threshold characteristic has been noted 
earlier by us [9] as well as by Fulton et al. [10]. No 
specific application was seen at the time. Intuition 
suggests, however, that this sinusoidal character 
should be reflected in the voltage developed across 
the SQUID, and calculations bear this out, as we 
shall presently show. 

The equations describing the threshold behavior 
of a dc SQUID are [9] 



Li 



Ii>C 



SQS^-1 



XI ; 



Figure 1— Diagram of a dc SQUID with bias current. The in- 
ductances and junctions on the two sides may be different. 
The crosses represent the junctions, including resistance and 
capacitance. 



cos <}>2=— cos <f>,/(a+j3 cos <f>i), (3) 

I m = /„, sin <J>, + 7 02 sin <f> 2 , (4) 

2tt<J>/<I>o=<J>2-(J>i+/32 sin (J> 2 -/3, sin <J>,, (5) 

where 4>i and 4> 2 are the quantum mechanical phase 
differences across the two Josephson junctions, 7 i 
and 7 02 are the respective critical currents, a=I 02 / 
7oi, 0i = 2irLiIoi/Q>o, P = 27rLI Q2 /<& , L =L { +L 2 , and 
L x and 7, 2 are the lumped inductances for the two 
segments of the SQUID. See figure 1. Equation (3) 
shows that there may not be a real solution for <\> 2 if 
the denominator on the right is less than unity. In 
fact if a and ft are each much less than unity, <f> 2 has 
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Figure 2 — The positive half of the 
threshold characteristic of a 
symmetric dc SQUID having 
L\ =1.2=2.0 pH and / ,=/ 02 = 
55/xA. 
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a solution only for <f>! very close to (2« + l)7r/2 
where n is any integer or zero. Choosing § x ~tt/2 
and noting that the requirements of a 4 1 and (i 4 
1 requires y3 2 4 1 (although /3] may still be of the 
order of unity), we find from eqs (4) and (5) that 



/ m ~/r 



,+/o2 sin (27rO/Oo-^2 sin <f> 2 + 7r/2+/3,) 

~/ OI +/o2COs(27ra>/0 +)3 I ). (6) 

This establishes the approximate sinusoidal behav- 
ior of the threshold characteristic of a sufficiently 
asymmetric dc SQUID with low inductance. 

In figure 3 we show the computed threshold be- 
havior of an asymmetric dc SQUID, having equal 
inductances L X =L 2 = 2.0 pH, and unequal critical 
currents 7 i = 100 juA, I 02 = 10 ft A. We also plot the 
perfect sinusoid of eq (6) for comparison. Calcula- 
tions show that the spectral purity of the threshold 
characteristic can be improved by making the in- 
ductances asymmetric as well. This remarkable pu- 
rity of the threshold characteristic is carried over 
into the voltage across the SQUID. 

The dynamic equations describing a dc SQUID 



are 






d 2 § 2 L x d§\ 
dt 2 ~ LlC2 dt 2 + R, dt 



Li d<j> 2 
R 2 dt 



+cf>i-<}>2+/3isin({>i-/3 2 sin<f> 2 + 277^-=0, (7) 



^,^1^,1^ 
Cl dt 2+C2 dt 2+ R l dt + R 2 dt 



+f 1 sin4> 1 +f^sin<k = 27r^. 

L,\ L, 2 Vq 



(8) 



Here the Cs are the junction capacitances and the 
R 's are resistances shunting the junctions. I b is the 
bias current injected as shown in figure 1. To cal- 
culate the voltage across the SQUID, we solve eqs 
(7) and (8) with a fourth-order Runge-Kutta tech- 
nique. We then time-average over typically 10 cy- 
cles to obtain the voltages shown in the figures. 
The precision of our calculations is estimated as a 
few parts in 10 5 . 

In figure 4 we show the calculated time-aver- 
aged voltage across the symmetric SQUID of fig- 
ure 2 for a bias current of 200 /xA, about twice the 
maximum critical current A sinusoid is also shown 
for visual comparison. Shunt resistances of 1 ft and 
junction capacitances of 0.66 pF are used for the 
two equal junctions. The junction capacitances are 
calculated by assuming a current density of 1000 
A/cm 2 , together with a specific capacitance of 12 
fiF/cm 2 , characteristic of Nb The voltage modula- 
tion at the bias current used is about 1 1 u,V relative 
to an average voltage of about 90 jxV, or about 12 
percent modulation. 

Figure 5 shows the frequency spectra in decibels 
for this case as well as for an asymmetric case. The 
dc component is suppressed. The spectra are ob- 
tained by calculating 32 voltage points in one pe- 
riod of the flux, then using a 32-point Fast Fourier 
Transform (FFT). Unity on the horizontal axis rep- 
resents the fundamental period <J> . The imprecision 
in the voltage calculations is amplified somewhat in 
passing through the FFT. Thus, points below 
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Figure 3— The positive portion of 
the threshold characteristic of 
an asymmetric dc SQUID with 
L,=L 2 = 2.0 pH, 7o, = 100 fiA, 
7 02 = 10 /iA. The solid line is the 
sinusoid of eq (6) with /3i=0.61. 
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Figure 4 — The time-averaged vol- 
tage across the symmetric 
SQUID of figure 2 with a cur- 
rent bias of 200 u,A (dots). 
Shunt resistances of 1 fl and 
junction capacitances of 0.66 pF 
are used. The solid line is a sinu- 
soid for visual comparison. 
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Figure 5 — Frequency spectra, in 
decibels, of the time-averaged 
voltage across two dc SQUIDs. 
The solid circles represent the 
symmetric SQUID of figures 2 
and 4. The crosses represent the 
asymmetric SQUID of figures 3 
and 6. The two SQUIDs have 
the same total inductance and 
critical current. Levels below 
about -70 dB are significantly 
affected by computational noise 
(see text). 



about —70 dB in figure 5 are significantly affected 
by computational noise. 

The solid circles in figure 5 show the spectrum 
of the symmetric SQUID of figure 4. The spectral 
purity of the voltage of this SQUID is better than 
that reported for an rf SQUID [4]. The second har- 
monic is about 26 dB below (5 percent of) the fun- 
damental, and the third harmonic is 50 dB down 
(0.3 percent). 

Figure 6 shows the time-averaged voltage vs. 
flux for the asymmetric SQUID of figure 3, again 
with a bias current of 200 uA and shunt resistances 
of 1 ft. The junction capacitances are calculated as 



above and have values of 1.2 and 0.12 pF, corre- 
sponding to critical currents of 100 /jA and 10 fxA. 
The calculations show an excellent voltage sinu- 
soid with a voltage modulation of 4 /xV relative to 
an average voltage of about 83 }xV, or about 5 per- 
cent modulation. A perfect sinusoid is also shown 
for visual comparison. The frequency spectrum of 
this SQUID is shown as the crosses in figure 5. 
The improvement in spectral purity of the asym- 
metric SQUID over the symmetric SQUID with 
the same total inductance and critical current is 
substantial. The second and third harmonics are 
now down 36 and 69 dB. The depth of modulation 
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Figure 6 — The time-averaged vol- 
tage across the asymmetric 
SQUID of figure 3 with a cur- 
rent bias of 200 uA (dots). 
Shunt resistances of 1 ft each 
are used. Junction capacitances 
are C\ = 1.2 pF and C 2 = 0.12 pF. 
The solid line is a sinusoid for 
visual comparison. 
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of the voltage is not as large as with the symmetric 
SQUID, however. We find this to be generally true 
— asymmetric SQUIDs have superior spectral pu- 
rity but less depth of modulation than symmetric 
SQUIDs with the same bias current and the same 
total inductance, critical current, and shunt resis- 
tance. 

The computed I-V curves of figures 7(a,b) show 
the differences in modulation for the two SQUIDs 
considered here. Curves are shown for O=0 and 



0=0.5<J>o. Note that these values of flux are not 
quite at the extrema of the voltage for the asym- 
metric SQUID, as figure 6 shows. Spectral purity 
generally increases as the bias current is increased, 
but at the expense of depth of modulation. An ex- 
ception occurs when an LC resonance significantly 
affects the current-voltage relation, as shown in 
figure 8. Here the shunt resistances are kept at 1 ft 
but the critical currents, and hence capacitances, of 
the SQUIDs considered above are doubled, and 
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Figure 7— Current-voltage characteristics for dc SQUIDs at 
<I>=0 (solid curves) and <J>=0.5<I>o (dashed curves). 
L\ =L 2 =2.0 pH and .R, =R 2 = 1 ft for both cases, (a) Symmet- 
ric SQUID with 7 01 =/ 02 =55 jxA, C, = C 2 =0.66 pF. (b) 
Asymmetric SQUID with 7 , = 100 /xA, 7 02 =10 fiA, C, = 1.2 
pF, C 2 =0.12pF. 
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Figure 8— Current-voltage characteristics for dc SQUIDs at 
*P=0 (solid curves) and *=0.5<I>o (dashed curves). 
L|=Z.2 = 1.0pHandR| = /?2=l ft for both cases, (a) Symmet- 
ric SQUID with 7 , = 7 02 = 110 /iA, C, = C 2 =1.32 pF. (b) 
Asymmetric SQUID with 7 i = 200 uA, / 02 = 20 u.A, C, = 2.4 
pF, C 2 =0.24 pF. 
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the inductances are halved to keep the LI C product 
unchanged. For the symmetric SQ UID — figure 
8(a)— an LC resonance at <V(277Vl,C,)=0.29 
mV is evident. (The resonance becomes sharper 
with larger values of shunt resistance.) Calculations 
for SQUIDs with apparent LC resonances show 
that the spectral purity of the voltage vs. flux 
curves is degraded. The asymmetric SQUID of fig- 
ure 8(b) does not show a prominent LC resonance 
because the two sides of the SQUID interfere; 
however, the larger capacitance of this example 
greatly reduces the modulation depth. It is thus im- 
portant that the junction capacitance be kept as 
low as possible so that the capacitive impedance 
does not shunt out the resistance. 

Perhaps the principal virtue of the low-induc- 
tance SQUIDs, and especially the asymmetric low- 
inductance SQUIDs, lies in the fact that the 
harmonics above the second are very greatly re- 
duced and are probably negligible. This is impor- 
tant because the earlier work with rf SQUIDs [5] 
found that the second harmonic could effectively 
be nulled, and that the higher harmonics consti- 
tuted the principal problem. That the second har- 
monic is also substantially reduced is of course an 
improvement. Since the basic source of systematic 
error in the measurement of rf attenuation is caused 
by harmonic distortion of the sinusoidal response 
of the system [4], the SQUIDs discussed here 
should have a distinct advantage over the rf 
SQUIDs. 

In 1982 M. Cromar of this laboratory made a 
preliminary study on the suitability of a thin-film 
dc SQUID as the detector element in an rf attenua- 
tor-calibrator system (unpublished). Using a resis- 
tively shunted symmetric SQUID whose LI C 
product was considerably larger than <£ , he 
showed that at low signal frequencies, the zeroes of 
the response approximated the zeroes of the J Bes- 
sel function closely enough that accuracy at rf fre- 
quencies comparable to that of the earlier 
measurements with rf SQUIDs might be realized. 
This research was not pushed further. 

The microfabricated dc SQUIDs are expected to 
have better thermal and temporal stability than the 
single-junction rf SQUIDs used earlier. The overall 
circuitry necessary to attain the output voltage 
should also be simpler with the dc SQUID. Imper- 
fections in microwave components were found to 
be a major contributor to the harn.onic distortion 
in the rf SQUID system [5]. Such components 
would not be used in the readout scheme for a dc 
SQUID. 

The 4-pH inductances used here are small, but 
SQUIDs with lower inductance have been fabri- 



cated. The low inductance of the proposed 
SQUIDs presents a potential problem of sufficient 
coupling, however. For example, if the mutual in- 
ductance between a 4-pH SQUID and the input 
line were only 4-pH, about 500 jllA would be re- 
quired in the line to produce one flux quantum in 
the SQUID. If one desires 200 nulls in the voltage 
response (100 flux quanta), about 50 mA must flow 
in the line in this case. Although higher than typi- 
cal, this value is still below what a superconducting 
stripline could support. Higher values of mutual in- 
ductance will decrease the maximum current 
needed. The mutual inductance to an input coil 
carrying the signal can in fact be made much larger 
than the self-inductance of the SQUID. Planar cou- 
pling to low-noise, low-inductance SQUIDs is the 
subject of much current research [11-13]. Of 
course, one may also choose to work with higher- 
inductance SQUIDs, accepting somewhat less 
spectral purity if the critical currents are kept at 
the same values, in order to decrease the maximum 
current needed in the stripline or to decrease the 
degree of coupling to the SQUID. 

The effect of noise is another consideration. It is 
beyond the scope of this note to undertake a de- 
tailed study of noise effects, which depend upon 
the readout method used. Ideally, of course, it is 
desirable that the noise of the SQUID plus its read- 
out system be dominated by the intrinsic noise of 
the SQUID. Readout schemes with this in view are 
under active investigation; Ketchen [11] discusses 
several of them. For a readout system operating at 
a frequency of about 100 kHz, \/f noise is not a 
consideration. However, even at lower frequencies 
where \/f noise dominates, a new readout scheme 
[14] for significantly reducing the noise has been 
developed. 

Properly fabricated dc SQUIDs have the lowest 
noise figures of any devices [11]. White noise de- 
creases with decreasing inductance, which is favor- 
able for the SQUIDs under consideration. The 
parameter 2tt kT/I <$>o is sometimes used to charac- 
terize the noise. If we take 7 =0.11 mA, which is 
the average of the critical currents used in the pre- 
ceding examples, we find the value 0.002 at 4 K for 
this noise parameter. The ultra-low-noise SQUIDs 
operate near this value. The spectral density S y of 
the voltage fluctuations in a resistor R is given by 
4kTR. For shunt resistances of about 1 ft, the 
voltage noise power spectrum is 2X 10~ 22 F 2 /Hz at 
4 K. The flux noise density S* is approximately 
S y /\dV/d<P\ 2 . The ultra-low-noise SQUIDs are 
built to be biased at that value of flux that gives the 
greatest energy sensitivity, which occurs at the 
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largest value of | aV/aO| . However, in the present 
case the flux will be swept over many periods of 
4> 0j so that each value of flux contributes almost 
equally. Since the voltage is nearly sinusoidal, 
| dV/a<l> | 2 varies almost as sine-squared, which we 
may replace by 1/2 times an amplitude, to a good 
approximation. Thus S^ — lkTR^l/iTrV^) 1 , where 
V is the voltage amplitude. For the symmet- 
ric SQUID of figure 4, F ^5.5 u-V so that 
5*~ I0~ n <$>l/Hz. These values are not far from val- 
ues quoted for low-noise SQUIDs [13-15]. Other 
expressions for the voltage noise power spectrum 
are available [14,15] which take into account noise 
contributions from the circulating current in the 
SQUID. The noise values from these expressions 
are not significantly different from those calculated 
above because of the very small inductance. Al- 
though these arguments do not establish that noise 
will not be a significant problem affecting the 
accuracy of determining the voltage nulls in rf 
measurements, they show that the SQUIDs pro- 
posed here share characteristics of low-noise 
SQUIDs, and thus are encouraging. 

The purpose of this note has been to suggest a 
new scheme for rf measurements using SQUIDs, 
and to demonstrate its potential advantages. A 
more extensive analysis of feasibility could include 
the following: simulations of the effect of noise 
upon the accuracy with which the zeroes of the 
response can be determined; studies of the magni- 
tude and effects of the parasitic inductance associ- 
ated with the shunt resistors; inclusion of the 
circuit carrying the signal — source, coupling coil, 
and load — to determine whether the nonlinear 
impedance reflected into the circuit by the SQUID 
is important; determination of the optimum choice 
of inductance, critical currents, asymmetry, and 
bias current to obtain the best combination of spec- 
tral purity and depth of modulation, or in short the 
greatest accuracy possible. 

The considerations presented here suggest that 
thin-film, low inductance, low critical current, dc 
SQUIDs, especially those with asymmetric critical 
currents, would be superior to rf SQUIDs in rf cur- 
rent, power, or attenuation measurements. 



This work has benefitted from discussions with 
several of my colleagues, to whom I express my 
thanks. 
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